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Abstract 

This  paper  presents  the  behaviour  of  a  proton  exchange  membrane  fuel  cell  (PEMFC)  connected  to  a  static  dc-dc  converter.  Two  different 
models  of  the  PEM  fuel  cell  are  obtained  from  static  measurements  and  impedance  spectroscopy.  The  paper  points  out  the  necessity  to  use  different 
models  according  to  the  type  of  study  performed  on  the  system.  The  comparison  of  models  at  high  semiconductors  switching  frequency  (25  kHz) 
is  illustrated.  Various  experimental  results  obtained  on  a  500  W  PEMFC  test  bench  are  compared  with  simulation  ones  to  illustrate  the  accuracy  of 
the  proposed  models. 
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1.  Introduction 

A  proton  exchange  membrane  fuel  cell  (PEMFC)  is  an  elec¬ 
trochemical  device  that  combines  hydrogen  fuel  and  oxygen  to 
produce  electricity,  heat  and  water  [1].  It  has  captured  world¬ 
wide  attention  as  a  clean  power  source  for  various  applications 
(electric  vehicles,  cogeneration,  . . .). 

In  many  applications  of  fuel  cells,  the  treatment  of  the 
delivered  energy  is  carried  out  by  static  converters.  The  inte¬ 
gration  of  fuel  cell  in  a  power  electronic  environment  (static 
converter-electric  load)  requires  the  knowledge  of  a  model  of 
the  generator.  This  model  allows  studying  the  whole  system  as 
presented  in  Fig.  1,  including  fuel  cell,  converter,  load  and  con¬ 
trol. 

Two  models  based  on  static  characteristic  and  electrical 
impedance  of  a  PEM  fuel  cell  are  presented  and  associated  with 
the  modeling  of  a  dc-dc  static  converter.  Then,  the  behaviour 
of  the  whole  system  is  considered,  first  for  low  dynamic  tran¬ 
sient  during  load  disturbance,  and  second  at  the  static  converter 
control  frequency  (25  kHz).  Different  experimental  results  are 
compared  with  simulation  ones  to  illustrate  the  proposed  mod¬ 
eling. 
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2.  Fuel  cell  model 


Two  different  models  are  considered.  The  first  one  is  obtained 
for  large  variation  of  the  load  parameters.  The  fuel  cell  is  rep¬ 
resented  by  means  of  its  voltage-current  characteristic  obtained 
in  static  operating  mode  [1].  The  second  ones  allows  studying 
fuel  cell  voltage  behaviour  at  the  semiconductor  switching.  This 
simple  model,  obtained  around  an  operating  point,  is  defined  by 
a  voltage  source  and  an  electrical  impedance. 

2.7.  Static  model  of  PEM  fuel  cell 


Fig.  2  presents  the  fuel  cell  voltage  as  a  function  of  the  cur¬ 
rent  magnitude  [1].  The  obtained  curve  is  composed  of  three 
main  regions  corresponding  to  the  electrochemical  activation 
phenomena  (region  1),  a  linear  part  where  the  voltage  drop  is 
due  to  electronic  and  ionic  internal  resistances  (region  2)  and 
the  last  part  where  the  diffusion  kinetics  of  gases  through  the 
electrodes  becomes  the  limiting  factor  (region  3).  This  last  zone 
is  characterized  by  a  sharp  voltage  fall. 

This  voltage-current  characteristic  fuel  cell  comportment  can 
be  defined  as  follow  [1]: 


Vpc  =  E  —  Aim 
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Fig.  4.  Transmission  line  model. 
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Fig.  5.  Proposed  model. 


where  Rd  is  the  diffusion  resistance,  Td  the  diffusion  constant 
time,  and  s  is  the  Laplace  operator. 

In  order  to  model  this  diffusion  convection  impedance  by 
discrete  elements,  it  can  be  first  considered  as  a  constants  dis¬ 
tributed  transmission  line  (Fig.  4). 

The  input  variables  are  linked  to  the  output  ones  by: 


Fig.  2.  Static  fuel  cell  characteristic. 

where  E  is  the  reversible  no  loss  voltage  of  the  fuel  cell;  /pc  the 
delivered  current;  iQ  the  exchange  current;  A  the  slope  of  the  Tafel 
line;  /l  the  limiting  current;  B  the  constant  in  the  mass  transfer 
term;  in  the  internal  current  and  Rm  represents  the  membrane 
and  contact  resistances. 

2.2.  Impedance  model  of  fuel  cell 
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with  m  =  y/rcs,  and  where  e  is  the  line  length. 

If  the  line  is  short  circuited  at  its  end  (x  =  e),  its  impedance 
becomes: 


Zi(s)  =  R 


In  the  case  of  a  PEM  fuel  cell,  the  electric  energy  is  obtained 
by  a  redox  electrochemical  reaction  [1].  The  simplest  small  sig¬ 
nal  model  for  this  type  of  reaction  proposed  by  Randles  can  be 
used  to  modelize  PEM  fuel  cell  (Fig.  3).  It  contains  a  transfer 
resistance  Rt,  which  characterizes  the  ions  transfer  phenomena 
at  the  electrodes  (only  cathode  is  represented),  a  resistance  Rm 
representing  the  membrane  and  contact  resistances,  a  double 
layer  capacitor  Cdi  and  a  diffusion  convection  impedance  Zws 
expressed  by  Nernst  for  a  finished  diffusion  length  by  the  fol¬ 
lowing  relation  [2,3]: 

th  (Jt& ~s) 
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where  R  and  C  are  the  total  line  resistance  and  capacitance  ( R  =  re 
and  C-ce). 

This  impedance  has  the  same  form  that  the  one  given  in 
relation  (2).  Then  one  can  represent  the  diffusion  convection 
impedance  by  such  a  model  with  a  finite  number  of  branches.  The 


PEMFC  stack 


Humidifier 


Fig.  3.  Randles  model. 


Fig.  6.  ZSW  (Ulm,  Germany)  500  W  PEMFC  test  bench. 
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Fig.  7.  Simplified  diagram  of  the  500  W  PEM  fuel  cell  system. 


simplest  model  proposed  is  represented  in  Fig.  5  and  contains 
for  this  impedance  two  different  branches,  R\-C\  and  R2-C2 . 

3.  Test  bench  presentation 

The  test  bench,  presented  in  Fig.  6  and  schemed  in  Fig.  7, 
uses  a  500  W  PEMFC  stack.  Air  is  supplied  from  a  compres¬ 
sor  through  a  humidification  unit,  and  pure  dry  hydrogen  comes 
from  bottles.  A  cooling  system  permits  to  impose  a  constant 
stack  temperature  by  heating  or  cooling  the  water  circuit.  A  pro¬ 
grammable  control  unit  contains  all  necessary  control  functions, 
such  as  reference  setting,  measurement  conditioning  and  secu¬ 
rity  shutdown.  The  stack  contains  23  cells  with  an  active  surface 
of  100  cm2  [6]. 

The  experimentation  is  carried  out  under  a  stack  tempera¬ 
ture  of  55  °C,  a  humidifier  temperature  of  45  °C,  and  a  use  of 
hydrogen  and  air,  respectively,  of  50  and  25%  at  atmospheric 
pressure. 

4.  Fuel  cell  model  identification 

4.1.  Static  characteristic  of  the  fuel  cell 

Experimental  identification  of  parameters  of  relation  (1) 
leads  to:  E=  -AG/2F  =  21.l  V  (for  T=55°C  and  23  cells), 
A=  1.35  V.  /'o  =  6.54  mA.  Rm  =  42  mQ.  B  =  1.19  V, /l  =  100  A  and 
in  =  230  mA.  Fig.  8  compares  the  experimental  stack  voltage 
with  the  calculated  ones  deduced  from  these  parameters. 

4.2.  Impedance  model  of  the  fuel  cell 

Thanks  to  the  impedance  spectroscopy  method  on  a  large 
range  of  frequency  [2-5]  and  a  least  square  method  identifica¬ 
tion,  the  parameters  of  the  model  of  Fig.  5  can  be  obtained. 

For  the  considered  PEMFC,  one  obtains  Rm  =  42m£2, 
R{  =  38  mQ,  Ccu  =  19  mF,  R\  =  40  mQ,  Ri  =  22  mQ,  C\  =0.15F, 
C2  =  2.37F.  Fig.  9  presents  the  measured  and  calculated 
impedances  by  the  proposed  model  (Fig.  5)  for  rated  current 
(40  A). 


Fig.  8.  Static  characteristic  of  the  500  W  PEMFC  (calculated  point  obtained 
from  relation  (1)). 


0.01- 

0 

-0.01  - 
-0.02- 
-0.03  - 
-0.04- 


Im  (Z) 


0.04 


0.08 


0.12 


0.16  Re(Z) 


» 

C* 


2  kHz 


\ 


0.05  Hz 


+  measured 
o  calculated 


Fig.  9.  Measured  and  calculated  impedance  for  the  rated  current  of  40  A  and  a 
frequency  range  of  0.05-2000  Hz. 


5.  Fuel  cell  converter 


5.1.  Converter  operation 

To  adapt  the  low  dc  voltage  delivered  by  the  fuel  cell  to  a 
standard  42  V  automotive  dc  bus,  a  boost  converter  -  operating 
in  continuous  mode  -  is  used,  as  shown  in  Fig.  10.  This  converter 
is  composed  of  a  high  frequency  inductor  L,  an  output  filtering 
capacitor  C,  a  diode  D\  and  a  main  switch  S\.  Switch  S2  is  a 


122 


I.  Sadli  et  al.  /  Journal  of  Power  Sources  156  (2006)  119-125 


Fuel  cell  Boost  converter  Load 


shutdown  device  to  prevent  the  fuel  cell  stack  from  short  circuit 
in  case  of  accidental  destruction  of  S\,  or  of  faulty  operation  of 
the  regulator. 


with: 
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When  Si  is  off,  state  differential  system  is: 

I  jX(t)  =  A2X(t)  +  U2(t) 

\  VBus (0  =  C2X(t) 

with: 


5.2.  Power  circuit  design 

The  design  of  the  converter  is  detailed  in  [7].  With  a 
12.5  V-40A  fuel  cell  and  assuming  an  efficiency  of  90%  for 
the  converter  with  a  switching  frequency  of  25  kHz,  a  bus  volt¬ 
age  ripple  of  2%  and  an  input  current  ripple  of  12%  of  the  rated 
value  one  obtains: 

Pout  =  450  W,  Ut  =  10.7  A,  /c,rms  =  19.1  A, 

P  —  72.9  [xH,  /s,rms  —  33.5  A,  /Dl,rms  =  21.9  A, 

C  =  357  pp  (5) 


Note  that  three  capacitors  C  (10  mF,  100  V,  7.32  A)  are 
selected,  associated  in  parallel  to  respect  the  constraint  on  RMS 
capacitor  current.  The  whole  input  capacitance  is  therefore 
30  mF,  so  that  voltage  ripple  is  much  lower  than  2%. 


5.3.  FC  converter  modeling 


5.3.1.  Switching  model 

In  order  to  design  controller  and  to  develop  system  simula¬ 
tion,  the  switching  model  of  fuel  cell  converter  has  to  be  known 
[8] .  The  non-linearity  due  to  the  power  semiconductor  devices  is 
negligible,  but  parasitic  resistances  of  the  power  semiconductor 
devices,  inductor  and  capacitor  are  considered.  The  load  of  the 
converter  is  classically  modeled  by  a  resistance  Rcu •  The  oper¬ 
ating  state  of  the  boost  converter  is  represented  by  the  a  state 
vector  X(t)  which  contains  the  fuel  cell  current  /fc(0  and  the 
capacitor  voltage  vc(t): 


he  (0 

vc(t) 


When  S\  is  on,  state  differential  system  is: 
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Rc  is  the  equivalent  series  resistance  (ESR)  of  the  output  capac¬ 
itor,  Rc  the  series  resistance  of  the  input  inductor,  Pm  the  whole 
drain-source  on  resistance  MOSFETS,  Vb  the  forward  voltage 
drop  of  diode  D\  and  bfc(0  is  the  fuel  cell  voltage. 


5.3.2.  Control  structure 

The  control  structure  of  the  fuel  cell  converter  designed  with 
a  classical  cascade  control  is  shown  in  Fig.  11. 

To  obtain  the  transfer  function  of  the  converter,  the  linearized 
differential  equations  from  relations  (6)  and  (7)  are  defined  as 
follow  [9,10]: 

d/pc(0  ~ 

L — - -  =  f’FcCO  -  Pl*’fc(0  -  (1  -  D)VBus(t)  +  Fbus^W 

J  at 

dfBus(t)  ~  ~  X’Bus(t) 

c =  (1  -  D)hc(t)  -  I¥C d(t)  -  -ppc. 

K  dt  RCu 

(8) 

where  D  is  the  nominal  duty  cycle  of  the  PWM  fuel  cell  con¬ 
verter,  d  the  duty  cycle  variation,  VbUs  the  nominal  dc  bus 
voltage,  fbus  the  dc  bus  voltage  variation,  he  the  nominal  fuel 
cell  current,  fpc  the  fuel  cell  current  variation,  Pch  the  load 
resistance  and  Pl  is  the  series  resistance  of  inductor  L. 

Note  that  series  resistance  of  capacitor  (Pc)  and  of  the  power 
semiconductor  devices  (Pm)  are  ignored.  From  Eq.  (8),  the  fuel 
cell  current  loop  can  be  modeled  by  the  following  open  loop 
transfer  function: 


I.  Sadli  et  al.  /  Journal  of  Power  Sources  156  (2006)  119-125 


123 


Fuel  Cell  Processor 


HEt 


H 


DC  Bus  Voltage 
Controller 


Air,  O 


Filter 


FCMea 


FC  current 
controller 


0 


1  FC 


< - i 


H:0 

— ► 


Heat 


Filter 


Fig.  11.  Cascade  voltage  and  average- current  control  structure  of  the  PEMFC  converter. 
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where  Vp  is  the  amplitude  of  the  PWM  saw  tooth  carrier  signal, 
K\  the  gain  of  measured  fuel  cell  current,  7fc  is  the  current  filter 
time  constant  and: 

^  _  (1  —  ^O^CH^FC  +  pBus 
7?l  +  ^ch(1  -  T>)2 
J  _  VbusRcrC 

(1  -  D)RcrIfC  +  ^Bus  ’ 


where  K 2  is  the  gain  of  the  measured  dc  bus  voltage,  7fv  is  the 
voltage  filter  time  constant  and: 


Gy  = 


Tzv  = 


h  = 
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(12) 


A  PI  controller  is  selected  for  dc  bus  voltage  controller,  and 
the  gain  Gcv  and  the  integration  time  constant  Tqv  are  set  to 
obtain  a  defined  phase  margin  (about  55°),  as  well. 

6.  Fuel  cell  with  static  converter 

6.1.  Load  variation 
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In  order  to  confirm  fuel  cell  converter  operation  and  to  com¬ 
pare  experimental  results  and  simulation  ones,  one  has  repre¬ 
sented  in  Figs.  12  and  13  the  system  response  to  a  load  step. 

(10)  The  model  used  for  the  PEMFC  is  the  static  ones  detailed  in 
relation  (1).  Fig.  13  shows  the  experiment  behaviour  of  the  fuel 


An  analogical  PID  controller  is  selected  for  fuel  cell  current 
controller.  The  derivation  time  constant  7qi  is  chosen  in  order  to 
compensate  the  filter  pole  -l/7fc.  By  classical  controller  design, 
the  gain  Gci  and  the  integration  time  constant  7q  are  set  to 
obtain  a  defined  phase  margin  (about  55°). 

As  far  as  the  dc  bus  voltage  loop  is  much  slower  than  the  fuel 
cell  current  loop,  the  open  loop  transfer  function  between  dc  bus 
voltage  measure  and  dc  bus  voltage  reference  can  be  written  as: 
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Fig.  12.  Simulation  PEMFC  converter  response  with  load  disturbance. 
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Fig.  13.  Experimental  PEMFC  converter  response  with  load  disturbance. 

cell  voltage  response  for  a  load  step,  this  response  is  almost  of 
resistive  nature.  The  voltage  drop  corresponds  to  low  frequency 
resistance,  which  is  identified  in  Fig.  9  at  5  mHz.  The  physical 
significance  of  this  resistance  is  due  to  the  sum  of  a  membrane, 
mass  transfer  resistance  and  diffusion  phenomena. 

(5.2.  High  frequency  behaviour 

Figs.  14  and  15  show  simulation  and  experiment  on  fuel  cell 
current  and  voltage  at  25  kHz,  the  switching  frequency  of  the 
boost  converter,  and  rated  power  (40  A,  12.5  V). 

The  simulation  and  experimental  results  follow  the  same 
behaviour.  The  static  model,  which  uses  a  dc  resistance,  leads 
for  a  given  duty  cycle  to  a  bigger  voltage  ripple  for  the  fuel  cell 
voltage.  The  ripples  obtained  with  the  dynamic  model  are  closed 
to  the  ones  deduced  from  the  experimental  results:  4.7  A  against 
4.2  A  for  the  fuel  cell  current,  and  211  mV  against  1 85  mV  for  the 
fuel  cell  voltage.  These  differences  are  not  so  important  if  one 
considers  that  one  has  taken  into  account  the  fuel  cell  model, 
and  the  parasitic  resistances  of  the  boost  converter.  The  rela¬ 
tionship  between  the  ripple  magnitude  of  fuel  cell  voltage  and 
current  gives  exactly  the  value  of  the  membrane  and  contact 
resistance.  This  value  is  identified  in  Fig.  9,  is  the  real  part  of  the 


high  frequency  (2  kHz)  measured  impedance,  and  is  estimated 
at  42  mQ. 

7.  Conclusion 
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Fig.  14.  Experimental  FC  voltage  and  current  at  40  A. 
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The  main  objective  of  this  work  is  to  use  different  mod¬ 
els  in  the  simulation  of  a  proton  exchange  membrane  fuel  cell 
with  dc-dc  converter.  The  experimental  and  simulated  results 
obtained  with  a  500  W  PEM  fuel  cell  are  closed  to  the  ones 
deduced  from  simulation  as  well  as  for  the  static  model  used 
in  low  dynamic  load  transient  than  for  the  dynamic  model  at 
the  converter  switching  frequency.  Experimental  results  at  high 
frequency  switching  of  the  semiconductor,  show  that  fuel  cell 
electrical  behaviour  is  mostly  resistive,  associated  to  membrane 
resistance,  which  is  accurately  described  by  the  dynamic  model. 
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